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O ' ABSTRACT 

\ We use high-resolution hydrodynamical simulations of large-scale structure formation to 

study the imprints of the loc al supe rclusters onto the full-sky Sunyaev-Zel'dovich (SZ) sig- 
nals. Following lMathis et alJ J2002I) . the initial conditions have been statistically constrained 
to reproduce the density field within a sphere of 1 10 Mpc around the Milky Way, as observed 
in the IRAS 1.2-Jy all-sky redshift survey. As a result, the positions and masses of prominent 
galaxy clusters and superclusters in our simulations coincide closely with their real counter- 
. parts in the local universe. We present the results of two different runs, one with adiabatic gas 

' physics only, and one also including cooling, star formation and feedback. By analysing the 

l/" - ) | full-sky maps for the thermal and kinetic SZ signals extracted from these simulations, we find 

C"^ ■ that for multipoles with I < 100 the power spectrum is dominated by the prominent local su- 

perclusters, and its amplitude at these scales is a factor of two higher than that obtained from 
unconstrained simulations; at lower multipoles (£ < 20) this factor can even reach one order 
of magnitude. We check the influence of the SZ effect from local superclusters on the cos- 
mic microwave background (CMB) power spectrum at small multipoles and find it negligible 
+3 , and with no signs of quadrupole-octopole alignment. However, performing simulations of the 

CMB radiation including the experimental noise at the frequencies which will be observed by 
the Planck satellite, we find results suggesting that an estimate of the SZ power spectrum at 
. large scales can be extracted. 
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1 INTRODUCTION 

It is well known that high-accuracy measurements of the 
anisotropies of the cosmic microwave background (CMB) are ex- 
tremely important for cosmology. By analysing their angular power 
spectrum, C(£), it is in fact possible to obtain tight constraints 
on the mechanisms from which the primordial fluctuations orig- 
inate and on the main cosmological parameters (or combinations 
of them) like the matter density, the baryonic density and the cos- 
mological constant contribution. For this reason, the astronomical 
community has been doing a big effort over the last years in order 
to improve the quality of the CMB data, both in terms of accuracy 
and resolution. 

The best available CMB measurements up to now were pro- 
vided by the Wilkinson Microwave Anisotropy Probe satellite 
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(WMAP; iBennett et al]|2003t) . The data obtained in its first year 
of activity gave strong support to the so-called concordance model, 
i.e. a spatially flat universe dominated by dark energy and cold dark 
matter, with nearly scale-invariant primordial fluctuations, almost 
gaussianly distributed. The WMAP result is confirming the general 
picture emerging by the analysis of other different datasets, like 
high-redshift supernovae la, weak lensing, galax y survey, L yman- 
q forest and galaxy clusters (see discussion in Sperge l et alj2003l 
and references therein). 

The statistical analysis and the theoretical interpretation of 
CMB data of so high quality necessarily require a good under- 
standing of all possible sources of contamination as well as the 
techniques to separate them from the CMB signal. The Sunyaev- 
Zel'dovich (SZ) effect iSunvaev & Zeldovichll972l . i.e. the Comp- 
ton scatter of the cold CMB photons by the hot electrons mainly 
located in the intracluster plasma, certainly represents one of the 
most important secondary anisotropies, i.e. temperature fluctua- 
tions generated along the light trajectory from the last scattering 
surface to the observer. Different studies have been devoted to the 
estimate of the SZ angular power spectrum which originates from 
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the whole galaxy cluster population (see, e.g., ICoorav et all2004 
and references therein). The results show that the SZ signal starts to 
dominate the cosmological one only on very small angular scales 
(I > 2000), while at the scales covered by WMAP is approxi- 
mately 3 orders of magnitude smaller. 

Even if its amplitude appears relatively small, the SZ sig- 
nal could be nevertheless significant at the sky positions cor- 
responding to nearby galaxy clusters. This has motivated the 
research of SZ imprints in the first-year WMAP temperature 
maps. This approach, mainly based on the cross-correlation of 
the CMB maps with the observed distribution of clusters and 
superclusters, as extracted in both optical and X-ray surveys, 
led to dif ferent claims for a positive detection at small angu- 
lar scales jM^ere - e^lJ|2004 |Hernandez^ontea gudo et alj|2004 
iFosalba & Gaztanagall2004 lAfshordi et alJl2004 . but only upper 
limits for the detectio n of t he SZ signal exist at la rger scales 
at present (Diego et al. 2003; Hernandez-Monteagudo et al. 2004; 
Iffirata et alJl2004 iHuffenberger et al J 12004 lHansen et aljEOoll) . 
In this sense the situation will largely improve in the near future 
thanks to the capabilities of the Planck satellite, whose launch is 
planned for 2007. In fact the range of frequencies covered by its 
receivers will be very large, and the noise level is expected to be 
very low. This will allow to extract the SZ signal for at le ast 10,000 
objects out to redshfit ?«1 (see, e.g. B artelm anrl200 lh . 

The CMB data from the WMAP satellite exhibit several 
anomalies at the lowest multipoles, among them the abnormally 
low v alue of the quadrupole amplitude (see, however, Efstathiou 
|2004 and the fact that the directions of t he quadrupole and oc- 
topole are unusually well alig ned (see, e.g.. lGaztanaga et all2003t 
de Oliveira-Cos ta et alj|2004 . There have been some speculation 
about the possibility that the SZ signal from structures in the local 
universe (the so-called local supercluster) could be the cause of the 
low quadrupole lAbramo & Sodreil2003h . Further, the quadrupole 
and oct opole directions happen to b e pointing toward the Virgo 
cluster Ide Olive ira-Costa et all2 004'). which opens the possibility 
of contamination from the local universe. 

The goal of this paper is to investigate the importance of the 
thermal and kinetic SZ effects produced by the local superclus- 
ters by comparing it with respect to the cosmological signal and 
to the contribution of more distant clusters. At this aim we will use 
the results of high-resolution hydrodynamical simulations realized 
starting from initial conditions l Math is et all2002|) which are con- 
strained to reproduce the density field recovered by the IRAS 1.2- 
Jy all-sky redshift survey. Using full-sky maps extracted from these 
numerical simulations, we will discuss the statistical properties of 
low multipoles in terms of amplitude and alignment. Moreover, tak- 
ing advantage of the different set of physical processes treated in 
the simulations we can assess the uncertainties related to the mod- 
elization. Finally we will perform simulations to discuss if the local 
SZ signal can be detected exploiting the capabilities of the Planck 
satellite. 

The plan of the paper is as follows. In Section 2 we describe 
the general characteristics of the constrained hydrodynamical sim- 
ulations of the local universe used in the following analysis. Section 
3 presents our method to produce full-sky maps of the thermal and 
kinetic SZ effects. Section 4 is devoted to the results of our anal- 
ysis. In particular we show the statistical properties of the maps in 
terms of angular power spectra and alignment of lower multipoles 
and we discuss the possibility of detecting the local SZ signal with 
future experiments, like the upcoming Planck satellite. Finally we 
summarize our findings and draw the main conclusions in the final 
Section 5. 



2 THE CONSTRAINED HYDRODYNAMICAL 
SIMULATIONS 

The results presented in this paper have been obtained by using the 
final output of two different cosmological hydrodynamical simula- 
tions of the local universe obtained starting from the same initial 
conditions but following a different set of physical processes. In 
more detail we used initial conditions similar to those adopted by 
iMathis et alJ 12002F) in their study (based on a pure N-body sim- 
ulation) of structure formation in the local universe. The galaxy 
distribution in the IRAS 1.2-Jy galaxy survey is first gaussianly 
smoothed on a scale of 7 Mpc and then linearly evolved back in 
time u p to z = 50 following the method proposed by Kolatt et al. 
1 1996) . The resulting fie ld is then used as a Gaussian constraint 
(Hoffman & Ribak 1991) for an otherwise random realization of a 
flat ACDM model, for which we assume a present matter density 
parameter Qo m = 0.3, a Hubble constant Ho = 70 km/s/Mpc 
and a r.m.s. density fluctuation as = 0.9. The volume that is con- 
strained by the observational data covers a sphere of radius ~ 110 
Mpc, centred on the Milky Way. This region is sampled with more 
than 50 million high-resolution dark matter particles and is em- 
bedded in a periodic box of ~ 343 Mpc on a side. The region 
outside the constrained volume is filled with nearly 7 million low- 
resolution dark matter particles, allowing a good coverage of long- 
range gravitational tidal forces. 

The statistical analysis made by Mat his et alJ 1200 21 demon- 
strated that the evolved state of these initial conditions provides a 
good match to the large-scale structure observed in the local uni- 
verse. Using semi-analytic models of galaxy formation built on top 
of merging history trees extracted from the dark matter distribu- 
tion in the simulation, they showed that the density and velocity 
maps obtained from synthetic mock galaxy catalogues have char- 
acteristics very similar to their observational counterparts. More- 
over many of the most prominent nearby galaxy clusters like Virgo, 
Coma, Pisces-Perseus and Hydra-Centaurus, can be identified di- 
rectly with haloes in the simulation, with a good agreement for sky 
positions and virial masses. In fact the positions of all identified 
objects differ less than the smoothing radius adopted in the initial 
conditions (7 Mpc), with the only exception of Centaurus, which 
is displaced by 9.6 Mpc. The agreement for the masses is within 
a factor of 2 which can be considered acceptable if one considers 
the uncertainties in the inferred observational masses, mainly when 
velocity dispersion is used. 

Unlike in the original simulation made bv lMathis etail 120021. 
where only the dark matter component is evolved, here we want to 
follow also the gas distribution. For this reason we extended the 
initial conditions by splitting the original high-resolution dark mat- 
ter particles into gas and dark matter particles having masses of 
0.48 x 10 9 Mq and 3.1 x 10 9 Mq, respectively; this corresponds 
to a cosmological baryon fraction of 13 per cent. The total num- 
ber of particles within the simulation is then slightly more than 108 
million and the most massive clusters will be resolved by almost 
one million particles. 

Our runs have been carried out with GADGET-2 ISpringell 

2005), a n ew version of the pa rallel Tree-SPH simulation code 
GADGET JSpringel et"al] 1200 ll). The code uses an en tropy- 
conserving formulation of SPH llSpringel & Hernquistll2002T) . and 
allows a treatment of radiative cooling, heating by a UV back- 
ground, and star formation and feedback processes. The latter is 
based on a sub-resolution model for the multiphase structure of 
the interstellar medium l Springel & Hernquist 2003). The code can 
also follow the pattern of metal production from the past history of 
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cosmic star formation iTornatore et al 12 0041. This is done by com- 
puting the contributions from both Type-II and Type-la supernovae 
and energy feedback and metals are released gradually in time, ac- 
cordingly to the appropriate lifetimes of the different stellar pop- 
ulations. This treatment also includes in a self-consistent way the 
dependence of the gas cooling on the local metallicity. 

As said, the results presented in this work are based on two 
different simulations starting from the same initial conditions. The 
first one, hereafter called gas, includes only non-radiative (i.e. 
adiabatic) hydrodynamics and has been originally analyzed by 
iDolae etall 120051) to study t h e prop agation of cosmic rays in the 
local universe. Hans en et alj «2005fl used the gas properties ex- 
tracted from this simulation to infer a prediction of the SZ effect 
from diffuse hot gas in the local universe. In that paper, the PSCz 
catalogue was used to map the matter density whereas the simula- 
tion was used only to assign a temperature for the gas with a given 
density. The SZ map realized in this way was cross-correlated to the 
WMAP data and upper limits were found. The second run, here- 
after called csf, exploits all the capabilities of the present version 
of the GADGET code, including cooling, star formation, feedback 
and metallicity; its outputs will be also used to study the detectabil- 
ity of the diffuse, warm intergalactic medium with future X-ray 
experiments (Kawahara et al. in preparation). 

The feedback scheme for the csf run assumes a Salpeter IMF 
and its parameters have been fixed to get a wind 
velocity of « 480 km/s. In a typical massive cluster the SNe (II 
and la) add to the ICM as feedback » 2keV per particle in an Hub- 
ble time (assuming a cosmological mixture of H and He); ~ 25 
per cent of this energy goes into winds. Note that these values can 
be considered an upper limit because a part of the ICM affected 
by the star processes could be at the present time out of the clus- 
ter virial radius. Moreover such feedback mechanisms cannot avoid 
the overcooling problem, commonly found in cosmological simu- 
lations: the corresponding overproduction of stars, acting mainly in 
the central regions, tends to amplify by a factor of » 2 the effects of 
the star population. Notice that the metal-dependence of the cooling 
function does not significantly change the global feedback prop- 
erties. Another signature that the feedback scheme is not enough 
efficient is the large amount of metals which is still locked inside 
the star particles: as a consequence, the resulting ICM metallicity 
is low, even if still compatible with observed values. A more de- 
tailed discussion of cluster properties and metal distribution within 
the ICM as resulting in simulations including the metal enrichment 
feedback scheme can be found in Tornatore et al. 120041) . 

The gravitational force resolution (i.e. the comoving softening 
length) of both simulations has been fixed to be 14 kpc (Plummer- 
equivalent), which is comparable to the interparticle separation 
reached by the SPH particles in the dense centres of our simulated 
galaxy clusters. 

In Table Q we report the main characteristics of the most 
prominent clusters in our simulations. In particular, we computed 
the (mass-weighted) temperature within one tenth of the virial ra- 
dius, and the central Compton-y parameter. In agreement with pre- 
vious results in the literature (see, e.g., Tornatore et al. 20q3), we 
find that the temperatures for a given object in the csf simulation 
are always larger than in the gas simulation by a factor of approx- 
imately 15-20 per cent. Including radiative cooling causes in fact 
a lack of pressure support, with a subsequent heating of the in- 
falling gas by adiabatic compression. The opposite trend is present 
for the central Compton-y parameter, which can be larger by a fac- 
tor up to 2 in the gas simulation than in the csf on e (see also 
IWhite et aj2002tlda Silva et al 120041: iMofl et aj2005T) . For com- 



pleteness, we also list the available observational estimates, which 
are in rough agreement. Notice that only for Coma we averaged 
in the simulation the value for Compton-j/ within an area of 450 
kpc to be compatible w ith the beam size of the observational data 
(Battistelli et al. 2003). Finally we would like to remark that the 
observed temperatures are not directly comparable with the simu- 
lation results, because they are extracted from X-ray spectroscopic 
data: as shown by Mazzotta et al ] fc004l) (see also lRasia et all2005t 
IVikhlininl2005h . there could be a significant bias produced by the 
complexity of the thermal structure. For completeness we added 
in the table the values for the spectroscopic-like temperature T SL , 
which approximates t he spectr oscopic temperature better than few 
per cent iMazzotta et alj2004l) . 



3 FULL SKY MAP MAKING 

Hot electrons, which are present in the intracluster medium, scat- 
ter by inverse Compton the cold photons of the cosmic mi- 
crowave background (CMB) radiation and re-distribute them to- 
wards higher frequencies. The result is the so-called SZ effect 
ISunvaev & Zeldovichl 19721) which originates a temperature decre- 
ment below 217 GHz, and an increment above. 

Given a direction 0, the change of the CMB temperature 
Tomb produced by the thermal SZ effect is given by 

ATcmb 



= y(0)g(x) , 



(1) 



TcMB 

where g(x) = x ■ coth(s/2) — 4 and x = hv/kTcMB- 

The Compton-y parameter is related to the three-dimensional 
thermal electron density, n e , and to the electron temperature, T, by 

ha 



y{d) 



dl n e (6>,Z) T(6,l) , 



(2) 



a T being the Thomson scattering cross section. 

The kinetic SZ effect is produced by the motion of the intra- 
cluster gas with respect to the CMB. The corresponding change of 
the CMB temperature can be written as 

ATcmb 



TCMB 

with 



-w(ff) 



w(6) = / dln e (8,l) v r (e,l) 



(3) 



(4) 



where v r is the radial component of the cluster velocity. 

The line of sight integrals which are present in the definitions 
of y (eq.|2j and w (eq.[4) can be computed in the simulations by ex- 
ploiting the SPH kernel. Here we employ the spline kernel defined 
as iMonaghan & Lattanzioll 985) 



W(r/h) = £ 



l-l(r/hf + l(r/h) 3 




< (r/h) < 1 

1 sC (r/h) < 2 ,(5) 

2 < (r/h) 



where h is the so-called SPH smoothing length and r represents the 
impact parameter. In this way eq.j2j becomes 



Vint 



E 

3 



Pj 



(6) 



where rrij and pj are the mass and density of the j-th gas particle, 
Uj is the corresponding Compton y-parameter and dj (r) is the pro- 
jected distance with respect to the position r along the line of sight. 
A similar relation can be written for Wmt- 



4 K. Dolag et al. 



Table 1. Comparison of the properties of the 6 most promin ent clusters, listed in C olumn 1 . The observed cluster temperatures T obs (Column 2) are taken from 
iMohr et alj |1999^. with the exception of the data for Virgo i Sanderson et al. 2003) and A3627 i Ikebe et al. 2002). Columns 3 and 5 refer to the mass-weighted 
temperature within 0.1 X R v i T , as computed in the gas and csf simulations (T^w and T c 1 ^ f w , respectively). Columns 4 and 6 refer to the spectroscopic-like 
temperature i Ma zzotta et all2004) within 0.1 X iJ v i r , as com puted in the gas said erf simulations (Tj^ and T^, respectively). In Column 7 we report the 
value for the Compton-?/ parameter, as obtained for Coma by Battistelli et al. 1 2003). Columns 8 and 9 refer to the values for y computed in the gas and 
csf simulations (y e zs and y ca i, respectiv ely). For Coma we averaged the Compton-j/ within an area of 450 kpc to be compatible with the beam size of the 
observational data i Battis telli et all2003t) . 

Cluster T obs [kev] T»[kev] Tj£[kev] T™[kev] T^\kev] y obs [10- 5 ] y e ^[W- 5 ] 2/ csf [10- 5 ] 



Coma 8.21+0.16 6A 58 73 73 9 - 6 ±i!i 92 64 

Virgo 2.55+°°g 3.5 3.0 4.1 3.5 4.8 3.1 

Centaums 3.54+ u ' § 3.8 3.6 4.6 4.5 6.7 4.2 

Hydra 3.10t ;" 3.2 3.2 4.5 4.6 9.3 4.2 

Perseus 6.33+°^ 5.8 5.2 6.1 5.8 9.5 6.8 

A3627 5.62^0^ 3.7 3.6 4.4 4.4 6.6 4.2 



In principle the sum has to be done over all particles, but in 
the case of a compact kernel where W(r /h) gets zero at large r/h 
(as the one we use) it can be restricted to those particles having a 
distance from the line of sight smaller than twice their smoothing 
length. 

In practice, however, one would like to obtain the value aver- 
aged over all lines of sight crossing the pixel, instead of the value 
along the line of sight crossing the centre of the pixel only. This 
can become a problem when the projected size of the structure pro- 
ducing the signal gets smaller than the pixel size. Therefore one 
usually applies the so called gather approximation, where, given a 
pixel, all particles - whose projections either overlap or completely 
fall inside the pixel - are taken into account. In this case the value 
of the Compton- y parameter can be written as the sum of the con- 
tributions of all relevant particles as: 

y pix = V ^VjNi^-W^idj/hi) , (7) 
^ Pi Aj 

3 

where dj is the projected distance between the particle and the pixel 
centre and Wmt is the integrated kernel. The area Aj associated to 
each gas particle can be approximated as the square of the cubic 
root of the corresponding volume, i.e. Aj = (mj/pj) 2 ' 3 . The nor- 
malization factor Nj is given by 

Y J NjA pix W int {dj/hj) = Aj (8) 

pix 

to ensure the conservation of the quantity y when distributed over 
more than one pixel. 

The method needs two further corrections. First, for pixels 
which are only partially overlapped by the projected SPH smooth- 
ing kernel, we approximate the associated area j4 p i x as 

A pix = ,/A~x (dj +0.5 yfA~- hj) . (9) 

To compute the contribution to such a pixel, the weight of the in- 
tegrated kernel Wint is then corrected by a factor A p i^/A p i^. This 
must be taken also into account when normalizing the integrated 
kernel. We checked that this correction leads to very small changes 
in the majority of the pixels, but its effect can become substantial in 
the high-signal regions corresponding to the cluster cores, specially 
when the particles are distributed over a small number of cells. 

Second, when a particle contributes to one pixel only, we fix 
the integrated kernel Wint to be unity. This is important when the 
area associated to such particle is much smaller than the pixel area 
and ensures that the value corresponding to the particle is com- 
pletely given to this pixel: in this case the normalization factor Nj 




Figure 1. Sketch of three different geometrical situations for particles con- 
tributing to one pixel (represented by the black square). Particle 1 (in red), 
having a smoothing length hi and being at a distance n from the pixel 
centre, falls completely inside the pixel, but does not overlap the centre. 
Particle 2 (in green) covers only partially the pixel. For particle 3 (in blue) 
the whole pixel lies within its projected radius. 

can become significantly smaller than unity to ensure the conserva- 
tion of the quantities. Finally we notice that, since the SZ effect is 
independent of distance, no distance factor has to be included into 
the previous equations. 

Figure Q exemplifies the three different situations previously 
discussed by showing how the gas particles can contribute to one 
pixel (here represented by the square). In such cases, the weight 
contributed by particle 1 (shown in red) is Ai/A p ^, the weight 
for particle 2 (in green) is W[ n t{r2 / 'h-£)A2 / A2 and the weight for 
particle 3 (in blue) is Wi n t{r2 /h2)A p i x /As. Note that for particles 
2 and 3 the weight needs to be normalised to unity when summed 
over all pixels to which the particles are contributing. 

Our maps have been constructed by applying the previous 
technique directly to the HEALPIx representation of the full sky 
(Gorski et al. 1998). We fix the nside parameter to be 1024, which 
results in an all-sky representation with » 1.2 x 10 7 pixels. In 
order to avoid spurious effects the line of sight integral has been 
performed between a minimum and maximum distance (r m i n and 
r max , respectively). In particular we fix r m i n — 5 Mpc to be sure 
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that the observer is outside the SPH smoothing radius of all parti- 
cles, and r max = 110 Mpc, which represents a conservative limit 
of our high-resolution region. 



4 RESULTS 

4.1 Full-sky maps 

In Fig. [2] we show the results of the application of the previous 
method to obtain a full-sky map of the Compton-y parameter in 
supergalactic coordinates. The upper panel corresponds to the sig- 
nal resulting from our gas simulation, covering the local universe 
up to 110 Mpc. It is possible to recognize the most prominent fea- 
tures: Perseus, Virgo, Centaurus, Hydra, A3627 and Coma (very 
close to the northern pole of the map). Also the presence of the fila- 
mentary structure connecting the largest clusters is still evident. In 
the lower panel we show the same map where we add the contri- 
bution coming from more distant objects. This has been computed 
by exploiting the results of Schaefer et al. (2004a,b), who build SZ 
signa l maps by using the Hubble volume simulation i Colbe rg et alJ 
2000; Jenkins et al. 2001) and suitably inserting the outputs of dif- 
ferent hydrodynamic re-simulations of single clusters. Notice that 
the physical processes considered in these re-simulations, i.e. non- 
radiative hydrodynamics, are the same included in our gas simu- 
lation, so no biases are introduced in combining the results. With 
respect to the Sc haefer et alJ l2004ll work we add only the objects 
having a distance from the observer larger than 110 Mpc, because 
the closer ones are already represented in our simulation. The re- 
sulting map appears more patchy, having a significant contribution 
coming from a large number of distant clusters. 

The corresponding maps for the kinetic SZ signal w are shown 
in Fig. [3] Again the most important structures of the local universe 
are evident in the upper map and continue to give the dominant 
contribution in the map including distant objects, shown in lower 
panel. 

It is worth to notice that within the simulation we measure a 
velocity at the observer' position of v x ~ —270 km/s, v y « 420 
km/s and v z ~ —350 km/s (assuming supergalactic coordinates). 
This leads to a absolute velocity of v ~ 610 km/s, which is in good 
agreement with the observed value (see, e.g., the WMAP analy- 
sis bv lBennett et all2003l) . Again, the observer' velocity points to- 
ward a direction which differs from the observed one less than 10 
degrees: this emphasizes once more that our simulations are giv- 
ing a good representation of the large-scale structure in our local 
universe. 

By analyzing the maps obtained from the results of the csf 
simulation (not shown here) we find only small differences with re- 
spect to the previous results, mainly in correspondence of the high- 
density regions, where the cfs simulation tends to give smaller sig- 
nals. This difference can be quantified by looking to the pixel dis- 
tribution as a function of y and w. This is shown in Fig|4| where 
the contribution of the local universe only is considered. Concern- 
ing the thermal SZ effect (left panel), the histograms are in general 
very similar, with only a slightly higher frequency for higher val- 
ues in the gas simulation. The differences are more evident in the 
pixel distributions for the kinetic parameter w. again high (abso- 
lute) values have higher probability in the gas simulation than in 
the c\s/one. 

Analyzing the gas bulk motions, we find no significant dif- 
ference in the gas and cs/runs. Moreover the winds driven by the 
feedback process are very efficiently stopped short after leaving the 
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Figure 5. The power spectra for the SZ effects. Comparison between the 
cosmological CMB power spectrum (upper black solid lined) and the SZ 
signals obtained in the different simulations: Hubble volume simulation 
(black lines), local constrained gas simulation (blue lines); local constrained 
gas simulation combined with the Hubble volume simulation for larger dis- 
tances (green lines). Upper and lower curves refer to the contributions from 
thermal and kinetic SZ effects, respectively. All spectra are taken in the 
Rayleigh-Jeans region where ATcmb/Tcmb = —2?/. 

star forming region and consequently are not contributing to the ki- 
netic SZ signal. Therefore the difference observed in Fig.|4]reflects 
the fact that the pressure in the csf simulation is smaller than in the 
gas one. Again the magnitude of the effect is in good agreement 
with previous findings 1 W hite et all20 02). 

Notice that the distribution for w is clearly non-gaussian, with 
a more extended tail toward negative values, produced by the clus- 
ter motion inside the local universe. Since the spectral signature of 
the kinetic SZ effect is the same of CMB, a perfect removal of this 
contribution could be difficult. 

It is also worth to notice that the mean value of the Compton-y 
parameter computed over the whole sky is dominated by the back- 
ground map extracted from the Hubble volume simulation. For the 
total map (Fig. [2] lower panel) we get a mean value of 3.3 x 10~ 7 , 
whereas we find a much smaller mean value (4.1 x 10 -8 ) when 
considering only the map from the constrained simulation (gas 
run). As discussed before, this value slightly decreases when us- 
ing the cfs simulation (3.2 x 10 -8 ). Note that all these values are 
well b elow the upper limit derived from COBE FIRAS data ( Fixsenl 
l200ll . 

4.2 The angular power spectrum 

In order to quantify the amount of signal produced by the SZ ef- 
fects, we compute from our maps the angular power spectrum 
C(£). In particular we estimate the spectra in the Rayleigh-Jeans 
region, where ATcmb/Tcmb = —2y. 

The results are shown in Fig|5| where we compare the power 
spectrum of the primordial CMB radiation (here represented by the 
upper black line corresponding to the best fit of the WMAP data 
obtained bv lSpereelet"aT] '2003) to the signals extracted from the 
numerical simulations here considered. In particular the plot shows 
C(£) for the gas simulation and distinguishes the contributions 
coming from the local universe from the one produced by more 
distant galaxy clusters. Both the thermal and kinetic SZ effects are 
considered. The results confirm that for £ < 2000 the cosmologi- 
cal signal is dominating the SZ one, being approximately 3-4 orders 
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Figure 2. Full-sky map of Compton-?/ parameter in galactic coordinates. The upper panel shows the results from the gas simulation, covering the local 
universe up to 110 Mpc from the Milky Way, while in the lower panel we added the contribution from more distant objects, as estimated by Schaefer et al. 
(2004a). The position of the most prominent structures is indicated in the maps. 
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Figure 3. The same as Fig[2| but for the kinetic SZ effect, shown by using asinh(u)) to better display positive and negative values. 
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of magnitude larger. Moreover the kinetic signal is always smaller 
by a factor of 10-20 than that produced by the thermal SZ effect. 
By comparing th e results of our simulation with the estimates ob- 
tained by Schae fer" et alj 120041), who populated the dark matter 
only Hubble volume simulation IColberg et all20olj|jenkins et alJ 
2001) with individual, adiabatic cluster simulations, we notice that 
the SZ effect from the local universe (objects with distance smaller 
than 110 Mpc from the Milky Way) is the most important contribu- 
tion up to scales corresponding to £ w 100. This is due to the larger 
angular size of the local supercluster structure. On the contrary, the 
signal at larger multipoles (£ > 200) seems to be dominated by 
the more distant sources, which are lying outside of our simulated 
volume. 

In order to quantify the contribution coming from the most 
prominent clusters, we compute C(£) from our maps by applying 
the galactic cut and masking the regions where the 7 largest lo- 
cal galaxy clusters are located. The resulting power spectrum (not 
shown here) changes drastically for I > 50 falling off much more 
rapidly, while remains unchanged for £ < 50. This is an expected 
result, because the richest clusters are contributing to the power 
around 2-3 degrees, whereas the total signal from all the rest is 
dominating the large-scale power. 

We further checked the contribution from diffuse gas by mask- 
ing out in the maps the pixels with Compton-y parameter smaller 
than a given value (y — 10~ 6 ): in such a way only the signal from 
cluster cores remains. The resulting power spectrum is not strongly 
affected, showing that the diffuse gas contributes minimally to the 
observed SZ power. In Fig.|7l we show the power spectrum includ- 
ing and excluding the diffuse gas contribution (solid and dotted 
lines, respectively) along with the power spectrum extracted from 
the Hubble volume simulation (dashed line). Note that the Hub- 
ble volume simulation gives a lower power spectrum at the largest 
scales when compared to the constrained simulation, even when the 
diffuse gas is removed from it. 

Having two different constrained realizations of the local uni- 
verse, the first one with adiabatic physics only, the second one in- 
cluding also cooling, star formation and feedback (the gas and csf 
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Figure 6. Effects on the power spectra of the inclusion of different phys- 
ical processes in the constrained simulations of the local universe. Blue 
and green curves show C; as extracted from the gas and csf simulations, 
respectively. Upper and lower panels refer to the thermal and kinetic SZ ef- 
fects, respectively. All spectra are taken in the Rayleigh- Jeans region where 
ATcmb/Tcmb = — 2y. 
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Figure 7. The large-scale power spectrum of the SZ effect obtained from the 
constrained simulations including (solid line) and excluding (dotted line) 
the diffuse gas contribution. The dashed line shows the power spectrum of 
the Hubble volume simulation. 



simulations, respectively), it is possible to discuss how the previ- 
ous results depend on the set of the included physical processes. 
This is done in Fig[6| where we compare directly the power spectra 
obtained from the two simulations. We find a higher signal in the 
adiabatic simulation, both for the thermal and kinetic SZ effects. 
The amplitude of this difference is very similar in both cases, being 
approximately a factor of 2. This is in agreement with the results 
obtain ed from individual cluster simulations (see, e.g. JWhite et alJ 
I2002I) . Therefore we do not expect the signal to get stronger when 
including additional physical processes and we can take the results 
obtained from the gas simulation as a robust upper bound. 



4.3 Alignment of lower multipoles 

Several anomalies have been dete cted at low m ultipo les in 
the WMAP data. In part icular iTegmark etai] J2003I) and 
Ide Oliveira-Costa et ail 120041) find the CMB quadrupole and oc- 
topole to have an unusually high degree of alignment, which is sig- 
nificant at the 2cr level. Moreover they show that the directions of 
the quadrupole and the octopole are close to the direction of the 
Virgo cluster which could suggest a possible link to structures in 
the local universe. By analysing our SZ templates, we find that 
the directions of the quadrupole and octopole (shown in Fig. [5] 
for the gas simulation) are very different from that extracted from 
the WMAP data: in our maps the quadrupole and octopole point 
toward (b,l) » (42.3°, -177.6°) and (b,l) « (32.5°, 88.1°), 
respectively, while the preferred axes in the WMAP data are in 
the directions (b,l) ~ (58.8°, —102.4°) for the quadrupole and 
(b,l) « (62.0°, -121.6°) for the octopole. Moreover the ampli- 
tudes of both quadrupole and octopole are far too small to have 
such an effect. 

Finally it has been claimed by different groups that the 
quadrupole o f the C MB is anomalously low [see the discussion 
in lEfstathioul f2004) and references therein] and that this could be 
cause d by structures present in the local universe lAbramo & Sodrd 
120031) . Again, we find that the amplitude of the SZ effect is far too 
low to have any relevant effect. 



4.4 Detectability of the local SZ signal with the Planck 
satellite 

We investigate now whether the local SZ effect can be observed us- 
ing the upcoming Planck satellite. The very low noise level and the 
huge frequency range of the Planck experiment could make it capa- 
ble of seeing this effect at large scales. Note that the CMB itself is 
the most important contamination when trying to look for the SZ ef- 
fect at large scales ( Hansen et all2 005). As the CMB has the same 
temperature at all frequencies, using a map produced by computing 
the difference between two frequency bands could be a promising 
way of detecting the SZ effect (this strategy was attempted for the 
WMAP data in lHansen et alJl20o3) . This makes the observation 
of the SZ effect limited only by the instrumental noise. Unfortu- 
nately galactic foregrounds also have several frequency dependent 
components, and these need to be eliminated at a very high degree 
in order not to confuse detections of the SZ signal. Here we will 
assume that these have been completely removed outside a conser- 
vative galactic cut. 

In order to estimate the sensitivity to the local SZ effect for the 
Planck experiment, we simulated Planck data for the 70, 100, 143 
and 217 GHz channels, by using the technical specifications given 
at the Planck home page 1 . In more detail, we apply the following 
procedure. 

• For each simulation, we generate a CMB sky and smooth it 
by using the beams corresponding to the 4 channels of interest. The 
smoothing is performed in the harmonic space by using beams with 
size of 14, 9.5, 7.1 and 5 arcmin for the channels at 70, 100, 143 
and 217 Ghz, respectively. As the band width for each frequency 
channel is expected to be much smaller than the frequency differ- 
ence between the channels, we approximate the band width to be 
infinitely narrow 

• We generate white, non-uniform noise for each channel and 
add it to the maps. The Planck experiment is expected to have 
highly correlated noise which would affect the noise power spec- 
trum at the largest scales. Since here we are only interested in a 
rough estimate of the sensitivity of the Planck experiment to the 
SZ effect, we neglect this effect. Consequently our results could be 
slightly over-optimistic at the largest scales. 

• We add the SZ template to each map by using the frequency 
dependence of the SZ effect. 

• We construct three different maps, corresponding to the dif- 
ference between the maps at 217 and 100 GHz (called map A), 
between 143 and 100 GHz (map B) and between 143 and 70 GHz 
(map C). 

• We apply the KpO galaxy and point source mask used by the 
WMAP team (and publicly available at the Lambda website 2 ) be- 
fore calculating the power spectrum of the difference maps. 

In Fig. [5] we show the results of a set of 300 different simu- 
lations (the number of simulations has been chosen to obtain error 
bars correct within a few percent). The solid line corresponds to the 
power spectrum of the template outside the KpO galactic cut. The 
three shaded areas refer to the la spread from simulations of the 
three difference maps, map A having the smallest spread, and map 
C the largest one. The huge difference between the different pairs 
comes from two different factors: first, the bigger the frequency 
difference, the bigger is the SZ effect in the resulting difference 



1 http://www.rssd. esa.int/index.php?project=PLANCK 

2 http://lambda.gsfc.nasa.gov/ 
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Figure 8. Comparison between the quadrupole (left panels) and octopole (right panels) obtained from the WMAP internal linear combination (ILC) map 
(upper panels) and from SZ template extracted from the gas simulation (lower panels). 



map; second, as the CMB is absent in the difference maps, the lim- 
iting factor is the noise level which is highly different in different 
frequency channels. We notice that assuming white noise and a per- 
fect foreground subtraction outside the KpO cut, we can expect the 
local SZ effect to be detected by the Planck satellite. Even if real 
error bars might be larger than our simulations show, one will still 
be able to compare data with the predicted local SZ effect. Finally 
note that analyzing difference maps is a useful tool for checking 
that foregrounds and other systematic effects of the experiment are 
well understood. For this purpose it will be of high importance to 
know well the expected local SZ effect as it will give an important 
contribution to the difference maps. 



5 CONCLUSIONS 

We used two hydrodynamical cosmological simulations which are 
designed to represent the observed large-scale structure of the local 
universe (up to approximately 110 Mpc from the Milky Way) to 
investigate the imprints of the extended SZ signal caused by local 
superclusters onto the cosmological full-sky signal. The two sim- 
ulations, which assume the standard ACDM model, started from 
the same initial conditions, but included the treatment of a different 
set of physical processes: the first one assumes adiabatic physics 
only, while the second one follows also cooling, star formation and 
supernovae feedback. 

We find that the largest and most prominent structures ob- 
served in the SZ map constructed from our constrained simulations 
are caused by local superclusters like Pisces-Perseus or the Cen- 
taurus supercluster region. Even Virgo, a relatively poor but nearby 
cluster, belongs to the most prominent and extended features of the 
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Figure 9. The estimated power spectra of three difference maps from 
Planck simulations. The solid line shows the power spectrum of the local 
SZ template. The shaded bands correspond to the la spread of the esti- 
mated power spectra in 300 simulations for the following difference maps: 
A = 217-100 GHz, B = 143-100 GHz and C = 143-70 GHz. Map A has 
the lowest spread and map C has the largest one. The KpO galactic cut has 
been used in these simulations. Note that all spectra are normalized to the 
expected SZ effect at 100 GHz. 



map. Comparing the power spectra of the thermal SZ signal with 
the on e derived from the Hubble volume simulation I Schaefer et al. 
2004), we find that for I < 100 these structures lead to an am- 
plitude which is at least twice as high (one order of magnitude for 
£ < 20). When comparing the derived power spectra for the kinetic 
SZ this factor is even higher. 
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However, the overall amplitude of the effect is far too small 
(3-4 orders of magnitude) to have any significant influence on the 
lowest multipoles of the CMB. There have been some claims that 
the anomalously low CMB quadrupole and octopole observed by 
the WMAP satellite could be explained by the imprints of local su- 
perclusters. We have shown that the amplitude of the SZ low mul- 
tipoles must be considerably much higher in order to produce such 
an effect as well as to have an influence on the observed quadrupole 
and octopole. Moreover, the directions on the sky of the quadrupole 
and octopole inferred from our simulations do not coincide with the 
observed direction in the WMAP maps. 

Both for our gas simulation as well as for the Hubble vol- 
ume simulation, radiative cooling processes and feedback by star 
formation have been neglected, which could change the absolute 
amplitude of the inferred SZ signal of galaxy clusters. However, 
using a simulation which includes all these effects (the csf simula- 
tion), we even find a slight decrease of the SZ signal, similar to the 
findings obtained by using isolated cluster simulations. This will 
make it even more unlikely to explain the abnormal lower multi- 
poles observed in the WMAP data with the SZ signal of the local 
universe. 

Nevertheless, we find that when taking advantage of the huge 
frequency range and low noise level of the Planck satellite, an 
estimate of the SZ power spectrum at large scales could be ob- 
tained. We have shown this by producing a set of simulated CMB 
maps with experimental noise at several frequencies. Analyzing the 
power spectrum of maps obtained by taking the difference of CMB 
maps at different frequencies, we show that the power spectrum of 
the SZ effect is limited only by experimental noise and foreground 
contamination. We have found that the predicted noise level of the 
Planck satellite is low enough to allow for the detection of the lo- 
cal SZ power spectrum at low multipoles provided that the galaxy 
can be removed with sufficient accuracy. It is out of the scope of 
this work to obtain an accurate estimate of the predicted detection 
level by the Planck satellite, for which correlated noise as well as 
an accurate study of foreground removal procedures is required. 
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